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By performing local projective measurements on a two-qubit entangled state one can certify in 
a device-independent way up to one bit of randomness. We show here that general measurements, 
defined by positive-operator-valued measures, can certify up to two bits of randomness, which is the 
optimal amount of randomness that can be certified from an entangled bit. General measurements 
thus provide an advantage over projective ones for device-independent randomness certification. 


The non-local correlations observed when measuring 
entangled quantum particles certify the presence of in¬ 
trinsic randomness in the measurement outputs in a way 
that is independent on the underlying physical realiza¬ 
tion of these correlations. While this relation between 
nonlocality and randomness had been noted by different 
authors since the seminal work by Bell mi, it is only 
recently that the tools to quantify the intrinsic random¬ 
ness produced in Bell setups were provided [2HU- These 
tools were initially introduced in the context of device¬ 
independent randomness generation m iMs], but have 
also allowed us to obtain a much better understanding of 
the relation between randomness and Bell violations, two 
of the most fundamental properties of quantum theory. 
For instance, today we know that maximal randomness 
can be certified from arbitrarily small amounts of non¬ 
locality or entanglement m, or that maximal random¬ 
ness certification is possible in quantum theory, but not 
in general theories restricted only by the no-signalling 
principle HU. 

Despite all this progress, there are still fundamen¬ 
tal questions on the relation between randomness, non¬ 
locality, and entanglement that remain completely unex¬ 
plored. In this work we consider and solve one of them: 
we obtain the maximal amount of randomness that can 
be certified in a standard Bell scenario involving local 
measurements on one entangled bit or ebit. In order 
to achieve this maximum the use of general measure¬ 
ment beyond projective ones, often known as Positive- 
Operator-Valued Measures (POVM), is necessary. Thus, 
our results and techniques are also interesting because 
they provide one of the few examples in the context of 
Bell non-locality where the use of these general measure¬ 
ments provides an advantage over standard projective 
measurements (other examples can be found in jl4lll5j ). 

We formulate the relation between randomness and 
non-locality in the setting of non-local guessing games 
as considered in (J. Such games consist of two users, 
Alice and Bob, and an adversary Eve. Alice and Bob 
perform local measurements on two separate quantum 
systems, labelled by A and B. There are nq and ms pos¬ 
sible measurements on particles A and B , each producing 


r a and possible results. Measurement choices are la¬ 
belled by x and y, with x = 1 ,.. ., mA and y = 1 , . .., ms , 
while the corresponding results are labelled by a and 
b, with a = 1 ,,r a and b = 1,..., re, respectively. 
The behavior of Alice and Bob’s systems is character¬ 
ized by the finite set of mA x ms x ta x tb probabilities 
P = {P(ab\xy)}, where P(ab\xy) is the probability that 
outcomes a and b are obtained when performing mea¬ 
surements x and y on particles A and B. In our non¬ 
local guessing game, P is assumed to be given, i.e. it 
is a promise on the behavior of Alice’s and Bob’s sys¬ 
tems. The aim is for Eve to guess as well as possible 
Alice’s outcome for a certain input x. To achieve this, 
Eve can prepare Alice’s and Bob’s system in any way 
compatible with the given behavior P and the laws of 
quantum physics. A strategy S for Eve consists in i) a 
tripartite quantum state |'I')ab£: on a composite Hilbert 
space PL a ®Hb®He of arbitrary dimension characteriz¬ 
ing the possible correlations between Alice’s, Bob’s and 
Eve’s system, ii) for each value of x, a POVM A x on Ha 
with elements A a \ x and for each value of y, a POVM B y 
on Hb with elements B b \ y characterizing the local mea¬ 
surements of Alice and Bob, Hi) a POVM Z on He with 
elements Z a whose result is Eve’s best guess on Alice’s 
outcome. Such a strategy is compatible with P if 

P(ab\xy) = ('P ABE\A a \ x 0 B b \ y 0 I\A> abe)- (1) 

The figure of merit of the game is the probability that 
Alice’s output and Eve’s guess coincide, maximized over 
all strategies S compatible with P (a set which we denote 
S P ): 

G(x, P) = max ABE\A a i £ 0 / 0 Z^^abe)- (2) 

oG«Sp z ' 
a 

We refer to this quantity as the local guessing probability. 

We can also introduce a variant of the game in which 
Eve attempts to guess both Alice’s and Bob’s outputs for 
a given pair of inputs (x, y) in which case her strategies 
involve a POVM Z with elements Z ab and the figure of 
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merit is 

G(x, y, P) = max 'V'^abe\A a \ s <g> B h \y <g> Z^I^abe)- 
oG«Sp z —' 
ab 

(3) 

We refer to this quantity as the global guessing probabil¬ 
ity. 

The local and global guessing probabilities quantify the 
predictability of the result of measurement x. or of pair 
of measurements x and y , by a quantum observer with 
an optimal description of the experiment. Taking minus 
the logarithm in base two of these quantities gives a mea¬ 
sure of randomness expressed in bits. A bound on these 
quantities is often a central element in the analysis of ac¬ 
tual randomness generation of expansion protocols, such 
as HE], where it directly determines (up to statistical 
corrections) the final amount of randomness generated. 

Note that we always have G(x, P) > max a P{a\x) and 
G(x,y,P) > max a b P(ab\xy) since a simple strategy is 
for Eve to simply guess the most probable outcomes of 
Alice’s and Bob’s measurements without exploiting any 
further detailed information about Alice’s and Bob’s sys¬ 
tems. However, in general a non-trivial strategy per¬ 
forms strictly better than these trivial bounds. Note 
that the guessing probabilities satisfy a convexity prop¬ 
erty in the sense that if P admits the convex decompo¬ 
sition P = Q\P\ hr term of quantum realizable be¬ 
haviors Pa, then G(x,P) > Y^\<l\G(x, P\) and similarly 
for G(x,y,P). This follows from the fact that a strat¬ 
egy that Eve can follow is to prepare Alice’s and Bob’s 
system with probability q\ according to the behavior P> 
and use the optimal guessing strategy associated to P\. 
In particular in the case in which the given correlations 
P can be described by a local model, i.e. can be writ¬ 
ten down as a convex mixture of deterministic behaviors, 
one has G{x) = 1 and G{x, y) = 1 for any measure¬ 
ment x and y. However, the violation of Bell inequalities 
does not necessarily imply that G(x) < 1 or G{x,y) < 1. 
Upper-bounds on the guessing probabilities can be com¬ 
puted using the Navascues-Pironio-Acin (NPA) hierarchy 
for quantum correlations )9j, as shown in 0151 - 

In this work, our goal is to compute the maximal ran¬ 
domness that one can certify from one ebit. That is, 
our goal is to identify the correlations minimising the 
guessing probability among all those attainable by mea¬ 
suring an entangled state equal in its Schmidt form to 
| (j> + ) = (|00) + 111))/V^2. The obtained quantity defines 
the optimal amount of randomness that can be certified 
in a device-independent way using an entangled qubit. 

Local measurements on such a state can always be 
viewed as POVMs acting on a qubit since the local 
Schmidt dimension is 2. In the case where the local 
qubit measurements are projective, it is known that 1 
bit of local randomness 0 and 2 bits of global random¬ 
ness mm can be certified from an ebit. This is also 
the maximum that can be achieved under such measure¬ 
ments, since a qubit projective measurement has only 
two possible outcomes. Beating those bounds thus re¬ 
quires considering more general measurements, beyond 


projective. 

We start by stating a rather straightforward observa¬ 
tion: no more than 2 log d bits of local randomness and 
4 log d bits of global randomness can be certified by mea¬ 
suring an entangled state of dimension d x d. This follows 
from the convexity property of the guessing probabilities 
mentioned above and the fact that a POVM acting on 
a space of dimension d can always be decomposed as a 
convex sum of POVMs of at most d 2 outputs which 
can evidently contain at most 2 log d bits of randomness. 
In the case of qubits, no more than 2 bits of local ran¬ 
domness and 4 bits of global randomness can be certified, 
i.e. twice as much than using projective measurements. 

Our main result is to construct two examples of qubit 
correlations saturating this bound on the local random¬ 
ness. They thus provide examples of optimal randomness 
certification from one ebit. In the first example we prove 
analytically that the local randomness is 2 bits, while in 
the second we have to resort to semidefinite programming 
(SDP) techniques. 

First optimal construction — Our first construction is 
based on non-local correlations obtained by measuring 
the two-qubit maximally entangled state \(j> + ) with mea¬ 
surements x = 1,2,3 on Alice’s side corresponding to 
cr x , <Ty, tr z and with measurements y = 1,..., 6 on Bob’s 
side corresponding to (a x ± a y )/\/2, ( <J X A <x z )/\/2, {er y ± 
a z )/\^2. These measurements are chosen so that they 
produce the maximal violation of the Clauser-Horne- 
Shimony-Holt (CHSH) inequality [T7| with all the pos¬ 
sible pairs of measurements on the first particle, that is 
B{ 1, 2; 1, 2) = B( 1, 3; 3,4) = B{ 2,3; 5, 6) = 2y/2, where 

B{i,j; k, l) = E ik + Eu + E jk - Ep 

and E ik = J2 a b(~l) ab P(ab\ik). Finally a four-output 
measurement y = 7 is included on the second particle. 
This seventh measurement by Bob plays a special role in 
our construction as it is the one used to certify the two 
random bits. We denote it by R , for random, and its 
measurement operators by R b , with b = 1,... ,4. This 
measurement is rather generic but has to satisfy two re¬ 
quirements: (i) it is extremal in the set of qubit measure¬ 
ments and (ii) its measurement operators are such that 
Tr (R b ) = 1/4, V6. An example of such a measurement is 
given by a POVM where R b = 1 /4| i/jb){4>b | and the Bloch 
vectors corresponding to \ip b ) point to the direction of a 
tetrahedron. 

The measurements described above define the behav¬ 
ior P which Eve has to reproduce (possibly using other 
quantum realizations of arbitrary dimension). Clearly, 
the four-output measurement R, when acting on half of 
a maximally entangled state, gives P(b\y = 7) = 1/4 for 
all b and thus G{y = 7, P) > 1/4. As mentioned previ¬ 
ously, this is only a trivial upper bound on the amount 
of intrinsic randomness, which is generally far from be¬ 
ing tight. However, we now prove that for the correla¬ 
tions P defined above the bound is tight and, therefore, 
G{y = 7, P) = 1/4 and hence two-bits of local random¬ 
ness can be certified from P. 
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To understand the main intuition behind the choice of 
state and measurements in our construction, the idea is 
to exploit the fact that, roughly speaking, the only quan¬ 
tum way of getting the maximal quantum violation of the 
CHSH inequality, also known as the Tsirelson’s bound, is 
by performing anti-commuting measurements on a two- 
qubit maximally entangled state. We refer to these quan¬ 
tum correlations, which are unique, as Tsirelson correla¬ 
tions. The correlations generated in the previous quan¬ 
tum setup contain three blocks of Tsirelson correlations 
for the different pairs of settings on particle A and cor¬ 
responding measurements on B. This suggests that the 
only state and measurements that could have produced 
these correlations should be, up to local unitary trans¬ 
formations, precisely Pauli measurements X, Y and Z 
on A acting on a two-qubit maximally entangled state. 
Now, these three measurements when acting on half of a 
maximally entangled state allow reconstructing any mea¬ 
surement implemented on the other half. In fact they 
remotely project particle B onto the eigenstates of these 
three observables, which are tomographically complete. 
Therefore, it should be possible from the observed cor¬ 
relations to reconstruct and certify the POVM elements 
implemented on B and conclude that they certify the de¬ 
sired amount of randomness. As we will see, this intuition 
is not entirely true due to a problem with complex con¬ 
jugation, but it is enough to prove the desired result. In 
fact, we believe this construction is interesting per se and 
may find applications in other problems of certification 
and self-testing of quantum devices. 

After providing this intuition, let us now show that 
the given setup indeed certifies two bits of local random¬ 
ness. As described earlier, any quantum strategy of Eve 
corresponds to a tripartite state ’)abe, a set of three 
observables Ai,A 2 ,A 3 for Alice, six two-output observ¬ 
ables Bi,... ,Bq and one four-output measurement B 7 . 
This strategy should reproduce the given correlations P, 
as expressed 0, and thus also the CHSH expectations 
B( 1, 2; 1, 2) = B( 1, 3; 3,4) = B{ 2,3; 5, 6) = 2/2. This im¬ 
plies, given the self-testing property of the CHSH inequal¬ 
ity m and following Mosca-McKague ED], that up to a 
local isometry I'S’abe) = |0 + ) ab\iP)A'B'e- In addition 
Ai\^abe) = (X A <g> Ia')\4> + )abW A'B'E and similarly 
A 3 | 'S’abe) = ( Z A ® Ia')\(I) + )abWa'B'e■ On the other 
hand, A 2 \'S>abe) = {Ya ® M A >)\<j) + ) ab\^)A'b>e where 
Ma' is hermitian and unitary. Basically this states that 
the three observables A\ 1 A 2i A 3 are necessarily the Pauli 
measurements A', Y, Z acting on Alice’s system, except 
for A 2 for which a correction Ma' is needed. This cor¬ 
rection reflects the fact that the optimal measurements 
leading to the three maximal CHSH expectations given 
above are only defined up to a complex conjugation (see 
[201J for a discussion of this point). 

Let us now determine the action of the POVM B 7 . 
For simplicity of notation, let us denote it R and the 
corresponding outcome operators Rb . The correlations 
between the outcomes of this POVM and Alice’s observ¬ 
ables should equal those of the ideal set-up defined ear¬ 


lier, as expressed in ([I]). This means that ('S’abe\ A ^ ® 
Rb ® II'S’abe) = (</ l ~| ( V ® P/,| <j >+), where A 0 denotes the 
identity operator. 

Let us now note that the ideal POVM elements Rb used 
in the definition of P can be written as Rb = Yl f j, r b cr i* 
where {cr M : \i = 0,1,2,3} is the basis of the four Pauli 
operators and rjf are complex coefficients defining the 
POVM. We then have (</> + |er M ® Rb\(/> + ) = r£, hence 
('k abe\A^ ® Rb® II'S’abe) = r%. 

On the other hand, since Rb acts jointly on systems 
BB ', without loss of generality we can decompose its 
operators as Rb = X/ 07 x®i?j, w h ere {c M : /z = 0,1,2,3} 
is the basis of the four Pauli operators on B and R? 
are arbitrary hermitian operators on B'. Inserting these 
expressions for R b and using the specific form of \S> A be) 
and A\,A 2 ,A 3 enforced by the CHSH constraints, we 
find (4>a'B'e\I ® R b ® I\4> A'B'e) = r% in the case n / 2 
and (iPa'B'e\M a > ® Rg ® I\ipA'B'E) = rg when /j = 2 . 

Introduce the normalized states \<p B f) = (Mj, ® I ® 
Z e )\il>A'B' e)/ yjq±,e, where M j, is the projector on the ± 
eigenspace of M A > and Z e is the projector corresponding 
to Eve’s outcome e (without generality we can assume 
Eve’s measurement Z to be projective). We can then 
write 

±,e 

= E [«+<e < ;+>e + <1 -,e (4) 

e 

for fj, = 0,1, 3 and 

r b =Y 9±,e(^ ±,e | M A' ® Pfc ® I\<p ±,e ) 

±,e 

= E i q +’ e f b’ + ' e - 9-,e< ; "’ e ] (5) 

e 

for fj. = 2, where we have defined the coefficients = 

(( / 3 ±,e |I ® R% ® I\(p ±,e ). Note that these coefficients de¬ 
fine a family of valid qubit POVMs R ±,e with operators 
R±’ e = X/ ^b ,±,e <V These POVMs simply correspond 
to preparing an ancilla system B' in the state \<p%f) and 
performing the POVM R on the joint system BB'. 

Now redefine fj}’ - ’ 6 as above but with the sign of f b ’~’ e 
changed. Then this also defines valid POVMs, which 
are just the the complex conjugates of R ±,e . With this 
redefinition, we can now write eqs. 0 and 0 as 

< = E^er/ ;± ’ e (6) 

±,e 

for fj, = 0,1,2, 3. We can interpret this as providing a 
convex decomposition for the ideal POVM R in term of 
the POVMs R ±,e with respective weights q±. e - But since 
this POVM is extremal, we must have rE’ e = r b for all 
±, e. In particular f°’ ,e = r b = 1/4 for all ±, e. 
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Finally, let us now rewrite the guessing probability for 
the input y = 7 with these notations. We have 

3 

G(y = 7, P) = ^(4/ abe\I <8> Rb ® Zp'S’ abe) 

6=0 

3 

= g> R b ® Z b \ i/jA'B'e) 

6=0 

= ^2q±,bf° b ’ ± ’ b = 1/4, 

6=0 

which provide the two announced random bits. 

Second optimal construction— Our second construc¬ 
tion to generate two random bits from a qubit is slightly 
simpler but the certification of randomness makes use of 
the numerical SDP techniques introduced in [3] based on 
the NPA hierarchy for quantum correlations j5]. 

The construction is based on the elegant Bell inequality 
introduced in |14j . It is defined in a scenario involving 
three measurements on Alice’s side and four on Bob’s. 
All measurements have two outputs and the inequality 
reads 

Pel = £-11 + £l 2 — E 13 — £ 14 + £21 — £22 + £23 — £24 
+ £31 ~ £32 — £33 + £34 < 6. (7) 

The maximal known quantum violation of the inequality 
is equal to 4-\/3 and is obtained with a maximally entan¬ 
gled state, projective measurements Ai = a x , A 2 = a y , 
A 3 = a z on Alice, while Bob’s projective measurements 
are defined by the four vectors of a tetrahedron. In Ap¬ 
pendix A, we show that this known quantum violation is 
in fact optimal, which gives a new Tsirelson-type bound 
for quantum correlations. 

We introduce a four outcome measurement R, but now 
on Alice’s side. As above, this measurement will be used 
to generate the two random bits. Given these configu¬ 
ration of measurements, we define the modified elegant 
Bell inequality 

4 

Pei = Pei - k ^2 p ( a = h b = +l\x = 4, y = i) < 6, (8) 

t=i 

where k is an arbitrary strictly positive constant. As the 
last term in the inequality is always negative, the bound 
follows from the bound on /3 e i - The same argument im¬ 
plies that the quantum violation cannot be larger than 
4-\/3. If we use the known optimal qubit settings, given 
above, the only way of getting this maximal violation is if 
the POVM elements of measurement R are anti-aligned 
with the four projective measurements on Bob’s side, so 
that all probabilities P(a = i,b = +l|a; = 4, y = i) are 
zero. But then, the corresponding measurement, when 
acting on half of a maximally entangled state, define two 
random bits on Alice’s side. The intuition, then, is that 
the maximal violation of the modified elegant Bell in¬ 
equality should certify the generation of two random bits 
for measurement R. 
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FIG. 1. Lower bound on the randomness (— log 2 of the guess¬ 
ing probability) as a function of the visibility v ranging from 
0.85 to 1.0. The correlations are of the form vq + (1 — u)r, 
where q are the quantum correlations yielding the maximum 
violation of the respective inequalities, and r denotes the com¬ 
pletely unbiased correlations that assign the same probability 
to each measurement outcome. The figure was obtained using 
level 2 of the NPA hierarchy and thus only represents a lower 
bound on the maximal randomness (note for instance that we 
do not recover the optimal values of 2 bits for v = 1). Better 
noise resistance than the one provided by these curves may 
thus be obtained by performing a more complex analysis. 


We used the numerical techniques in 13] to bound 
the randomness present in the correlations maximally 
violating Recall that these techniques are based 

on SDP and, therefore, one has control over the pre¬ 
cision of the numerical result. Using these techniques 
at order 2 + AAB + ABB with an arbitrary-precision 
solver (23| 03 ] , we can show that the generated random¬ 
ness by measurement R in the previous setup is larger 
than 1.99999989474702 bits. 

Noise robustness — The numerical approach of [3| also 
allows one to study the robustness of the previous con¬ 
structions against noise. A typical noise model consists of 
mixing, with weights v and 1— v, with 0 < v < 1, the ideal 
quantum correlations with uncorrelated noise in which 
all outputs have the same probability. By decreasing v , 
often known as visibility, the amount of certified random¬ 
ness decreases. In Figure [I] we plot a lower bound on the 
generated randomness as a function of the visibility for 
the two previous constructions. It can be seen that the 
gain provided by the POVM is fragile, in the sense that 
a small fraction of noise, of the order of 0.01, makes the 
obtained randomness smaller than one bit, which is the 
randomness provided by projective measurements. These 
considerations are relevant, for instance, when thinking 
of a possible experiment showing the advantage of using 
POVM’s for randomness certification. A natural open 
question opened by our work is thus to identify robust 
setups for randomness generation using POVM’s. 

Global randomness — Before concluding, we would like 
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to briefly discuss the problem of global randomness. The 
question is whether it is possible to find Bell setups in¬ 
volving a pair of maximally entangled qubits allowing the 
generation of 4 bits of randomness. We mainly leave this 
question for future work. Nevertheless, we made some 
preliminary numerical searches, using slightly more com¬ 
plex variations of the previous construction based on the 
elegant Bell inequality. These constructions are described 
in Appendix B. As shown there, they can be used to cer¬ 
tify the generation of 2.8997 bits, which is both higher 
than the global randomness that can be certified with 
projective measurements and the local randomness that 
can be certified with general measurements. 

Conclusions — We have shown that an ebit can certify 
the presence of more than one bit of randomness locally 
and more than two bits globally. For the case of local 
randomness, we have found two constructions that can 
certify 2 bits of randomness, the maximal possible value. 
Both constructions involve a maximally entangled state 
and the three Pauli measurements on Alice’s side. In 
the first construction, the violation of three CHSH in¬ 
equalities are used to self-test the maximally entangled 
state and these three Pauli measurements on Alice’s side. 
This in turn allows to self-test a four-outcome extremal 
POVM on Bob’s, which generates the two bits of local 
randomness. Our second construction is instead based 
on a single inequality - the elegant inequality introduced 
in [14] , In this case, we had to resort to SDP techniques 
to put a bound on the local randomness. A possible way 
to proof this bound analytically, would be to first proof 
that the elegant Bell inequality also provides a self-test 
for the maximally entangled state and the three Pauli 
measurements on Alice’s side. It would then be possible 
to add the four-outcome measurement on Bob’s side and 
use the same proof as above to conclude that it generates 
two random bits. 

Both constructions are quite sensitive to noise, as a 
fraction of noise of the order of 0.01 makes the obtained 
randomness smaller than one bit, which can already be 
obtained with more robust constructions based on pro¬ 
jective measurements. A natural open question is thus to 
identify robust and optimal setups for randomness certi¬ 
fication using POVM’s. 

Finally, we also found a construction based on POVMs 
which yields more than two bits of global randomness, the 
best that can be obtained with projective measurements, 
but less than the theoretical maximum of four bits. It 
remains an open question whether this maximum can 
actually be attained. 
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Appendix A: SOS decomposition of the elegant Bell 
inequality 

The aim of this Appendix is to prove a tight upper 
bound on the quantum violation of the elegant Bell in¬ 
equality /3 e i defined by ([T]) in the main text. 

Using the maximally entangled two-qubit state \<fi + ) = 
(|00) + 111))/a/ 2 and projective measurements A\ = cr x , 
A 2 = a y , A 3 = a z on Alice’s side, while the following 
four projective measurements on Bob’s side 

Bi = ( a x - (T y + a z )/\/S 

B 2 = (cr x + dy - cr z )/Vs 

B 3 = (— -o y - cr z )/V3 

B 4 = ( -a x + a y + a z )/V?>i (Al) 

we obtain the Bell violation 4-\/3. 

We now show that this amount of Bell violation for 
inequality f3 e 1 is optimal. To this end, let us redefine the 
Bell operator, p e i = 4-\/31 — p e i, were /3 e i is the stan¬ 
dard Bell operator defined by the inequality ([7]) and the 
measurement operators by Alice and Bob: 

Pei = Ai (g> (Hi + B 2 — B 3 — R 4 ) 

+ A 2 <S> (Bi — B 2 + B 3 — B 4 ) 

+ A 3 <S> (Hi — B 2 — H 3 + H 4 ). (A2) 

Finding the largest quantum violation of the elegant Bell 

inequality ([7]) corresponds to maximising (minimising) 
the largest (smallest) eigenvalue of $ e i (p e i) over all quan¬ 
tum observables Ai and Bj by Alice and Bob. 

Using the fact that (A.j)tA.j = (BjpBj = 1, it is easy 
to see that p e i admits a sum of squares (SOS) decompo¬ 
sition in the form /3 e 1 = (\/3/2) ]T) A H A H\, where P\ are 
linear combinations of the operators Ai and U ; as follows 

Hi = (Ai + A 2 + A 3 )/y/3 - Hi 

P 2 = (Ai — A 2 — A 3 )/V3 — B 2 

P 3 = ( — Ai + A 2 — A 3 )/y/?> — H 3 

H 4 = (— Ai — A 2 + A 3 )/y/?> — H 4 , (A3) 

This implies that p e \ is positive semidefinite, hence 
Pei = 4-^1 — p e \ < 4-^/3!. Therefore, the maximum 
quantum value of P e 1 is upper bounded by 4^/3. Since 
this value is attained with a quantum realization pre¬ 
sented above, we have a new Tsirelson-type bound for 
quantum correlations. Note that by applying the meth¬ 
ods of Refs. [ 21 , 22 ] one gets the same upper bound of 
4-\/3 on the maximum quantum value of the Bell inequal¬ 
ity Pei. 
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Appendix B: Certifying global randomness 


Here we give a detailed description of the Bell setup 
used to generate global randomness exceeding 2 bits from 
two-qubit entangled states. Note that according to the 
main text it is not possible to certify more than 4 bits of 
global randomness from any two-qubit state. Since the 
observed statistics provide an upper bound to the intrin¬ 
sic randomness of the setup, in order to obtain 4 bits, we 
have to find a pair of 4-outcome extremal measurements 
on a two-qubit state which yield P(a,b) = 1/16 for all 
results. To this end, let us pick the maximally entan¬ 
gled state | (j) + ) = (|00) + 111))/a /2 and Alice’s and Bob’s 
rank-1 extremal POVM’s as follows: 


Ra =(1/4)(1 + u a ■ (?) 

Rb =(1/4)(1 + v b ■ a), (Bl) 


where the four outputs are a, b = {1,2, 3,4}, and tf de¬ 
notes the vector of three Pauli matrices (cr x ,a y ,a z ). 
Alice’s Bloch vectors u a are distributed as 


u\ = c( — 1, 5, 0) 

U 2 = c( —1, —8, 0) 

u 3 = c( 1, 0, -S) 

u 4 = c( 1, 0, 5), 


(B2) 


while Bob’s Bloch vectors Vb look as follows 


where 0 < 5 < 1. It is known that for a two-setting 
Bell inequality £n + 8 E 2 1 + £12 — SE 22 the quantum 
maximum is 2-s/l + S 2 . Using this result, it is straight¬ 
forward to establish the upper bound of 3 + Sx/1 + S 2 
on the maximum quantum value of the above inequality 
/9fork- This value can in fact be attained by using the state 
|</+) and observables A -4 = B± = a x , A 2 = — B 2 = er y , 
A 3 = B 3 = a z , and A x+3 = —v x ■ a for x = 1,..., 4 and 


B y+ 3 = — u y ■ tf for y = 1,... ,4, where the Bl och vec - 
tors above are defined according to equations ( B2|B3 1. 
In this construction, the diagonal term En + E 22 + E 33 
forces Alice and Bob to choose correlated measurements, 
whereas the other parts of the inequality are used to steer 
Alice and Bob’ s measu rements to the ones defined by the 
Bloch vectors ( B2|B3 ). 

We now introduce a four-outcome measurement on Al¬ 
ice’s and Bob’s side, which will be used to generate the 
glo bal r andomness. With the help of the above inequal¬ 
ity (B4), we define a modified Bell inequality similarly to 
equation <[sj) in the main text as follows 


4 

/3fork = /3fork - k^2,P(a = i,b = +l\x = 8 ,y = i + 3) 
2=1 


4 

-kJ2 P (a 

2=1 


+1 ,6 = = z + 3, y = 8 ) < 11, 

(B5) 


vi = c( — 6 , 0 , — 1 ) 

v 2 = c( <5, 0 , - 1 ) 

v 3 = c( 0 , - 6 , 1 ) 

V 4 = c( 0, <5, 1) 


(B3) 


where c = 1 / \J\ + S 1 2 3 is a normalization factor and S is 
a small positive constant. Note that the measurements 
defined above fulfill two necessary conditions to approx¬ 
imate the value of P(a, 6 ) = 1/16 for all results, (i) 
Vectors u a and Vb become orthogonal when <5 tends to 
zero, (ii) Both sets of vectors u a and Vb span the three- 
dimensional space (i.e. the POVM elements are complex 
valued). 

Let us now construct a 7-setting 2-outcome Bell in¬ 
equality which reads 


where A: is a strictly positive constant. The subtracted 
terms above will take care of al igning Alice and Bob’s 
four-outcome measurements (Bl) in the correct way (i.e. 
force them to take a form having Bloch vectors (B2) and 
(B3)). Hence, in principle the probability P(a,b ) = 1/16 
can be approximated for all a, b pairs provided 5 tends 
to zero. 

Solving the above problem at order 1 + AB of the NPA 
hierarchy we get a 121-dimensional moment matrix. By 
taking S = 0.0676946 and k = 1 , the problem in this level 
can be solved on a desktop, and our technique certifies 
2.8997 bits of global randomness. 


/3fork = -Eii + E 22 + E 33 
+ E\4 + 8E 2 4 + £15 — 8E 23 
— Eiq + 8E 3 q — En — 5 E 3 7 
+ 8E41 + £43 — 5 E 5 i + £53 
+ 8Eq 2 — £g3 -|— 8E^2 — £73 < 11 , (B 4 ) 
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